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ABSTRACT 

The  next  system  design  under  consideration  for  the  Los 
Alamos  Foil  Implosion  Project  (Trailmaster)  is  projected  to 
deliver  in  excess  of  15  MA  of  electrical  current  produced  by 
high-explosive  driven  flux  compression  generators  to  a  foil 
load.  A  plasma  flow  switch  is  being  investigated  as  the  final 
pulse  shaping  step  in  this  system.  The  performance  of  these 
switches  is  being  evaluated  using  a  wide  variety  of  computa¬ 
tional  tools  including  zero-,  one-  and  two-dimensional  MHD 
codes  and  a  3-D  view-factor  radiation  transport  code.  We 
are  concerned  with  the  effects  of  radiation  from  the  switch 
plasma  prior  to  switching  current  to  the  load,  and  the  inter¬ 
action  of  the  switch  plasma  on  existing  perturbations  and 
as  a  source  of  perturbation  on  the  imploding  load. 

I.  Introduction 

Researchers  at  Air  Force  Weapons  Laboratory  (AFWL) 
have  reported  considerable  success  in  efficiently  transferring 
current  to  an  imploding  liner  using  a  plasma  flow  switch.1 
In  the  Los  Alamos  foil  implosion  project,  Trailmaster,  we 
are  interested  in  this  approach  for  several  reasons.  Our  use 
of  explosively  driven  magnetic  flux  compression  generators, 
while  providing  an  inexpensive  source  of  high  current,  has 
led  to  a  considerable  pulse  shaping  problem.  The  plasma 
flow  switch  has  proven  more  efficient  in  transferring  cur¬ 
rent  than  other  closing  switch  options.  In  addition,  the 
plasma  flow  switch  protects  the  load  region  from  high  volt¬ 
ages  generated  by  an  opening  switch  until  the  current  is 
large  enough  to  provide  magnetic  insulation. 

In  the  Trailmaster  project  we  anticipate  using  plasma 
flow  switches  at  higher  current  levels  than  have  been 
switched  heretofore.  We  are,  therefore,  using  a  variety  of 
computational  tools  to  evaluate  switch  designs  for  the  Trail- 
master  project.  We  are  conducting  simulations  with  a  0-D, 
horizontal,  slug  model  to  find  the  optimum  plasma  mass 
and  flow  channel  parameters.  Our  1-D,  Lagrangian,  code, 
RAVEN,  will  run  in  planar  geometry.  It  is  being  used  to 
provide  a  variety  of  benchmarking  calculations,  radiation 
flux  for  ablation  studies,  and  initial  conditions  for  our  2-D, 
Eulerian,  calculations.  Our  3-D  view  factor  and  radiation 
transport  code,  GOPHER,  is  being  used,  in  conjunction 
with  RAVEN,  to  examine  the  effects  of  radiation  from  the 
flowing  plasma.  Finally,  we  are  in  the  process  of  making 
2-D,  Eulerian,  calculations  to  evaluate  the  effects  of  the 
flowing  switch  plasma  on  the  quality  of  the  liner  implosion. 


Zero-Dimensional  Modeline 


Our  0-D,  horizontal,  slug  model  is  included  in  the  elec¬ 
trical  circuit  package  in  the  RAVEN  code.  We  have  bench- 
marked  the  model  by  comparing  to  1-D,  planar  geometry, 
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calculations  with  RAVEN.  As  an  example  of  our  use  of  this 
model,  the  1.5  MJ,  Pegasus  capacitor  bank2  is  embarking 
on  a  series  of  plasma  flow  switch  experiments  to  assist  in 
designing  such  a  switch  for  Trailmaster.  We  have  used  the 
0-D  model  to  determine  the  optimal  mass  for  this  study 
such  that  the  switching  plasma  will  be  over  the  load  region, 
6.5  cm  down  the  channel,  at  the  peak  of  the  bank’s  quarter 
cycle.  This  quarter  cycle  period  is,  of  course,  a  function  of 
the  time  dependent  inductance  of  the  switch  itself.  Figure 
1  shows  the  0-D  model’s  predictions  for  current  and  plasma 
position  for  a  35  mg  plasma.  This  would  be  approximately 
30%  of  the  mass  in  the  switches  fielded  by  AFWL. 

III.  One-Dimensional  MHD  and  Radiation  Ablation 
Calculations 

We  are  using  RAVEN,  1-D,  in  planar  geometry,  to 
benchmark  the  0-D  model,  provide  radiation  flux  estimates, 
and  initial  conditions  for  2-D,  Eulerian,  calculations.  The 
radiation  induced  ablation  calculations  start  with  the  cal¬ 
culated  flux  output  from  RAVEN  and  view  factor  values 
from  the  GOPHER  code.  These  values  are  convolved  and 
the  resultant  flux  is  imposed  on  the  outer  most  zone  of  an 
exposed  surface  in  a  RAVEN  calculation.  In  making  these 
ablation  studies  we  were  concerned  that  the  interaction  of 
the  aluminum  plasma  in  the  switch  and  the  plastic  barrier 
foil  might  lead  to  higher  temperatures  than  would  be  pre¬ 
dicted  by  a  pure  aluminum  calculation.  We  have,  therefore, 
carried  out  comparison  calculations  with  pure  aluminum 
and  aluminum  plus  a  barrier  foil. 
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Fig.  1.  Slug  model  prediction  for  the  position  of  the  switch 
plasma  (dash-dot  line,  left  axis),  and  current  from 
the  Pegasus  capacitor  bank  for  a  plasma  mass  of 
35  mg. 


69  9 


Amps 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

JUN  1989 


2.  REPORT  TYPE 

N/A 


3.  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

Calculational  Evaluation  Of  Plasma  Flow  Switches  For  The  Los  Alamos 
Foil  Implosion  Project 

6.  AUTHOR(S) 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Los  Alamos  National  Laboratory  P.O.  Box  1663  Los  Alamos,  New  Mexico 
87545 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S )  AND  ADDRESS(ES )  10.  SPONSOR/MONITOR' S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002371.  2013  IEEE  Pulsed  Power  Conference,  Digest  of  Technical  Papers  1976-2013,  and 
Abstracts  of  the  2013  IEEE  International  Conference  on  Plasma  Science.  Held  in  San  Francisco,  CA  on 
16-21  June  2013.  U.S.  Government  or  Federal  Purpose  Rights  License. 

14.  ABSTRACT 

The  next  system  design  under  consideration  for  the  Los  Alamos  Foil  Implosion  Project  (Trailmaster)  is 
projected  to  deliver  in  excess  of  15  MA  of  electrical  current  produced  by  high-explosive  driven  flux 
compression  generators  to  a  foil  load.  A  plasma  flow  switch  is  being  investigated  as  the  final  pulse  shaping 
step  in  this  system.  The  performance  of  these  switches  is  being  evaluated  using  a  wide  variety  of 
computational  tools  including  zero-,  one-  and  two-dimensional  MHD  codes  and  a  3-D  view-factor  radiation 
transport  code.  We  are  concerned  with  the  effects  of  radiation  from  the  switch  plasma  prior  to  switching 
current  to  the  load,  and  the  interaction  of  the  switch  plasma  on  existing  perturbations  and  as  a  source  of 
perturbation  on  the  imploding  load. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OF  PAGES 

SAR 

4 

19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


In  these  comparison  calculations  we  have  used  Pega¬ 
sus  capacitor  bank  parameters  but  AFWL  mass  values.  In 
running  RAVEN  in  planar  geometry  we  approximate  the 
cylindrical  shape  of  the  channel  as  a  rectangle  with  width 
equal  to  the  circumference  of  the  center  of  the  channel  and 
height  equal  to  the  separation  of  the  electrodes.  In  this  cal¬ 
culation  we  have  used  a  polyurethane  equation-of-state  from 
the  SESAME  tables  and  not  allowed  it  to  become  conduc¬ 
tive.  The  wire  array  is  approximated  by  an  aluminum  foil 
with  the  correct  mass. 

The  positions  of  the  Lagrangian  zones  as  functions  of 
time  for  both  calculations  are  shown  in  Fig.  2.  Note  that 
in  both  calculations  the  upstream  (back)  zones  of  the  alu¬ 
minum  (in  the  direction  of  the  capacitor  bank  and  at  higher 
values  of  Z)  blow  off  but  are  caught  by  the  J  X  B  force  and 
collide  back  with  the  main  body  of  the  plasma  in  the  1.5  to 
2  n s  time  frame.  The  exact  time  is  dependent  on  the  zon¬ 
ing  resolution  which  is  finer  on  the  back  of  the  plasma  in 
the  aluminum  only  calculation.  The  collision  between  the 
aluminum  plasma  and  barrier  foil  starts  at  1.1  /is  but  its 
effects  are  not  felt  on  the  back  of  the  plasma  until  well  af¬ 
ter  the  2  /is  point.  The  temperatures  of  the  zones  in  both 
calculations  are  shown  in  Fig.  3.  The  highest  values  are  for 
those  zones  on  the  back  of  the  plasma,  where  the  current  is 
being  carried.  The  actual  values  are  quite  similar  between 
the  two  calculations.  The  barrier  foil  does,  however,  clearly 
reduce  the  spread  in  the  total  plasma.  In  comparing  the 
results  shown  in  Figs.  2  and  3  we  have  concluded  that  for 
most  purposes  calculating  with  aluminum  only  (no  barrier 
foil)  for  the  switch  plasma  is  sufficient. 
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Fig.  2.  Calculated  positions  as  functions  of  time  for  the 
Lagrangian  zones  with  (top)  and  without  (bot¬ 
tom)  a  separate,  polyurethane,  barrier  foil.  Total 
mass  is  the  same  in  each  calculation,  absolute  po¬ 
sitions  have  no  significance. 


0  0-3  1  L5  2  2-5  i  J-5  4 


Time  OxSec) 

Fig.  3.  Calculated  temperatures  for  the  Lagrangian  zones 
with  (top)  and  without  (bottom)  a  separate,  poly¬ 
urethane,  barrier  foil. 

One  possible  future  system  in  the  Trailmaster  project 
would  combine  a  Mark  IX  generator,3  an  explosively  formed 
fuse  in  a  flux-conserving  geometry,4  and  a  plasma  flow  switch. 
As  a  preliminary  test  for  such  a  system,  a  dynamic  load, 
similar  to  a  plasma  flow  switch,  has  been  designed.  A  blue¬ 
print  of  this  load  is  shown  in  Fig.  4. 

We  have  used  the  RAVEN  code  to  estimate  the  per¬ 
formance  of  such  a  system.  It  predicts  that  the  current  in 
the  dynamic  load  will  reach  approximately  16  MA  in  2.3  /is 
if  the  mass  of  the  plasma  is  the  same  as  that  used  in  the 
AFWL  experiments.  The  highest  temperatures  are  on  the 
“back”  of  the  plasma,  the  generator  side.  This  tempera¬ 
ture  reaches  8  eV  during  the  initial  shock  when  the  J  x  B 
force  overcomes  the  expansion  of  the  plasma  and  drives  the 
back  zones  of  the  plasma  into  the  main  body  of  the  plasma. 
During  most  of  the  calculation  this  temperature  is  closer  to 
6  eV.  We  are  concerned  that  the  radiation  from  this  6-8  eV 
plasma  might  lead  to  radiation  induced  ablative  channel 
closure. 

We  have  used  the  3-D  view  factor  and  radiation  trans¬ 
port  code  GOPHER,  to  model  the  flow  of  radiation  up  the 
channel  past  the  first  radiation  baffle.  We  have  concen¬ 
trated  our  attention  on  the  surface  indicated  by  asterisks 
in  Fig.  4.  At  this  point  the  channel  is  still  narrow  but  this 
surface  does  have  some  direct  line-of-sight  with  the  initial 
position  of  the  wire  array.  We  convolve  the  result  of  the 
view  factor  calculation  with  the  predicted  radiation  flux 
from  the  RAVEN  calculation  and  impose  the  resultant,  time 
dependent,  flux  on  this  surface  in  another  RAVEN,  planar 
geometry,  calculation.  The  predicted  flux  from  the  back  of 
the  load  reaches  a  peak  of  2.5  xlO12  W/m2  during  the  tem¬ 
perature  peak  at  341.08  /is.  During  most  of  the  calculation 


700 


IV.  Two-Dimensional  MHD  Calculations 


Fig.  4.  Propose  the  dynamic  load  for  future  Trailmaster 
system.  Surface  modeled  by  the  ablation  calcula¬ 
tion  is  marked  by  asterisks. 

the  flux  is  approximately  1  x  1012  W/m2.  We  performed 
calculations  assuming  that  the  channel  was  clean,  pure 
aluminum,  and  another  with  the  aluminum  coated  with  a 
0.1  mm  layer  of  a  hydrocarbon.  Far  more  ablation  occurred 
in  the  hydrocarbon  coated  calculation  and  the  time  depen¬ 
dent  positions  of  the  zone  boundaries  are  shown  in  Fig.  5. 
The  ablated  plasma  moves  about  0.8  cm,  less  than  halfway 
across  the  channel.  In  addition,  this  is,  almost  certainly,  a 
worst  case  calculation.  In  making  our  view  factor  calcula¬ 
tion  we  used  the  initial  position  of  the  wire  array  through¬ 
out.  By  the  end  of  the  4  ps  shown  in  this  calculation  the 
wire  array /plasma  would  have  moved  more  than  6  cm  down 
the  channel.  We  conclude  that  radiation  induced  ablation  is 
not  a  serious  threat  to  cause  channel  closure  in  this  design. 


Once  specific  designs  have  been  determined  for  our 
work  with  the  Pegasus  capacitor  bank  and  explosive  gen¬ 
erator  systems,  the  1-D  calculations  will  be  used  to  estab¬ 
lish  detailed  initial  conditions  for  the  2-D,  Eulerian,  MHD 
calculations.  Prior  to  this  we  are  using  such  2-D  calcula¬ 
tions  with  generic  initial  conditions  to  examine  the  interac¬ 
tion  between  the  switch  plasma  and  the  imploding  load.  We 
are  particularly  concerned  that  this  interaction  may  impose 
perturbations  on  the  load  that  could  initiate  a  magnetically 
driven  Rayleigh-Taylor  instability.  We  note  that  the  radia¬ 
tion  pulses  from  the  AFWL  experiments  have  a  FWHM  of 
the  order  of  0.2  /is  and  full  base  width  exceeding  0.7  ps.1 

Figure  6  shows  results  of  a  2-D  calculation.  In  this  cal¬ 
culation  we  are  using  a  constant  current  drive  of  10  MA. 
The  initial  radius  of  the  load  foil  is  4  cm.  In  our  2-D  work 
we  start  with  the  load  somewhat  expanded  and  in  a  plasma 


Fig.  5.  One-dimensional,  Lagrangian  calculation  of  radi¬ 
ation  induced  ablation  from  the  surface  indicated 
in  Fig.  4.  The  material  being  ablated  is  from  a 
thin  hydrocarbon  layer. 


Fig.  6.  This  figure  shows  the  interaction  of  the  switch  plasma  and  the  implod¬ 
ing  load  at  three  times  as  calculated  by  a  2-D,  Eulerian  code.  The 
switch  plasma  appears  to  initiate  a  Rayleigh-Taylor  instability  with  a 
wavelength  that  is  twice  the  height  of  the  load. 
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state.  In  this  case  the  density  of  the  load  aluminum  is  5.5  x 
10-s  gm/cm3  at  a  temperature  of  3  eV.  The  figure  shows 
density  profiles  in  the  region  of  the  interaction  between  the 
switch  and  the  load  at  the  indicated  times. 

The  results  shown  in  Fig.  6  clearly  indicate  that  the 
interaction  between  switch  plasma  and  the  load  can  ini¬ 
tialize  an  instability  in  the  imploding  load.  This  instabil¬ 
ity  starts  with  a  wavelength  that  is  twice  the  height  of  the 
load.  Therefore,  this  effect  may  be  exacerbated  by  the  cate¬ 
nary  which  tends  to  be  present  in  our  ultrathin  foils  and  is 
not  present  in  this  calculation.  These  results  also  suggest, 
however,  that  it  may  be  possible  to  ameliorate  these  prob¬ 
lems  by  tailoring  the  load  foils  so  that  there  is  more  mass 
near  the  region  of  first  contact. 
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